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SUMMARY 


A two-step hybrid analysis technique is presented for predicting the nonlinear 
steady-state temperature distribution in structures and solids. The technique is 
based on the successive application of the regular perturbation expansion and the 
classical Bubnov-Galerkin approximation. The functions associated with the various- 
order terms in the perturbation expansion of the temperature are first obtained by 
using the regular perturbation method. These functions are selected as coordinate 
functions (or temperature modes) and the classical Bubnov-Galerkin technique is then 
used to compute their amplitudes. 

The potential of the proposed hybrid technique for the solution of nonlinear 
thermal problems is discussed. The effectiveness of this technique is demonstrated 
by means of three numerical examples, which include the effects of conduction, con- 
vection, and radiation modes of heat transfer. 

Results of the study indicate that the hybrid technique overcomes two major 
drawbacks of the classical techniques: (1) the requirement of using a small parame- 

ter in the regular perturbation method and (2) the arbitrariness in the choice of the 
coordinate functions in the Bubnov-Galerkin technique. Therefore, the proposed tech- 
nique extends the range of applicability of the regular perturbation method and 
enhances the effectiveness of the Bubnov-Galerkin technique. 


INTRODUCTION 

In recent years, increasing attention has been devoted to the application of 
approximate analytical techniques to nonlinear heat-transfer problems. Analytical 
techniques have the major advantage over numerical discretization techniques of pro- 
viding physical insight into the nature of the solution of the problem. Moreover, 
analytical techniques can be used in conjunction with a partitioning scheme for the 
thermal analysis of individual components of practical (complex) structures. Two of 
the commonly used approximate analytical techniques are the regular perturbation 
method and the Bubnov-Galerkin technique. Review of the many applications of these 
techniques to thermal problems is given in a survey paper (ref. 1) and a monograph 
(ref. 2) . 

The regular perturbation method consists of the development of the solution in 
terms of unknown functions with preassigned coefficients. The unknown functions are 
obtained by solving a recursive set of differential equations which are, in general, 
simpler than the original governing differential equation of the problem. By con- 
trast, in the Bubnov-Galerkin technique, the temperature is sought in the form of a 
series of a priori chosen coordinate functions (or temperature modes) with unknown 
coefficients . 

Despite their usefulness in solving nonlinear thermal problems, the aforemen- 
tioned two techniques have a number of drawbacks. Regular perturbation techniques 
have two major drawbacks. The first stems from the fact that as the number of terms 
in the perturbation series increases, the mathematical complexity of the differential 
equations builds up rapidly. Therefore, for practical applications, the perturbation 
series has to be restricted to a few terms. The second drawback is the requirement 



of restricting the perturbation parameter to small values in order to obtain solu- 
tions of acceptable accuracy. The main shortcoming of the Bubnov-Galerkin technique, 
from a practical point of view, is the difficulty of selecting good coordinate func- 
tions (or modes). 

The aforementioned drawbacks of the regular perturbation technique have been 
recognized and a number of remedial actions were proposed. These included the use of 
a small number of terms (e,g,, two or three) in the perturbation expansion and either 
(1) generating "mimic functions" which give accurate numerical estimates of the solu- 
tion over the entire physical domain (ref, 3), or (2) applying a nonlinear transfor- 
mation (e.g,, Shanks transformation, ref, 4) to estimate the solution as the number 
of terms goes to infinity. However, as shown in reference 5, the success of these 
methods cannot be guaranteed in general, and the remedial actions may fail to produce 
satisfactory results. 

The successful experience obtained with the hybrid numerical procedure which 
combines the finite-element method with the Galerkin technique in reference 6 raises 
the question as to whether the use of a hybrid technique which combines both the 
standard regular perturbation method and the classical Bubnov-Galerkin technique 
might overcome the major drawbacks of the two techniques and provide a more effective 
approximate analysis procedure. The present study focuses on this question. Specif- 
ically, the objectives of this paper are (1) to present a hybrid perturbation/Bubnov- 
Galerkin technique for nonlinear thermal analysis and (2) to demonstrate the effec- 
tiveness of the proposed technique by means of numerical examples. The technique is 
described as a formal procedure without any attempt to justify it rigorously. How- 
ever, the numerical results presented herein are intended to give the analysts some 
insight into the potential of the proposed technique and to stimulate research and 
development of the mathematical foundations necessary to realize this potential. 

To sharpen the focus of the study, discussion is limited to nonlinear steady- 
state thermal analysis with continuous temperature fields in the space domain. Non- 
linear conduction, convection, and radiation modes of heat transfer are considered. 
The hybrid technique presented herein is also expected to be particularly useful for 
nonlinear transient thermal problems. 


SYMBOLS 

A cross-sectional area of fin 

B boundary differential operator (see eq, (2)) 

c circumference of fin cross section 

C,C^,C2f^2 constants defined in appendix 
e error norm defined in equation (17) 

h convective heat-transf er coefficient 

fe thermal conductivity coefficient 

thermal conductivity coefficients at T = T^ and T = 0 
L fin length (also, side length of plate) 
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differential operators defined in equations (1), (3), and (4) 

differential operators defined in equations (5), (6), (9), (10), (11), 
and (12) 

number of points at which numerical values of temperature are evaluated 
total number of coordinate functions 
range of indices (see eq, (8)) 
perturbation parameters 

functions defined in equations (10) and (12) 
temperature 

coordinate functions (temperature modes) defined in equations (7) 
and (8) 

amplitude of sinusoidal temperature variation (see fig. 1) 

Cartesian coordinates 
parameters (see table I) 
boundary of structure (or solid) 

nonlinear conductivity coefficients (see table I) 

constant (see table I) 

emissivity of. fin surface 

nondimens ional temperature (see table I) 

first eigenfunction for zeroth-order perturbation equation 
dimensionless coordinate functions defined in appendix 
radiation-conduction parameter (see table I) 
dimensionless coordinates (see table I) 

= 1-5 

Stef an-Boltzmann constant 
amplitude of eigenfunction 

undetermined coefficients (amplitudes of coordinate functions) 
domain of structure (or solid) 
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Subscripts: 


a ambient 

b fin base 

s effective sink 

BASIC IDEA OF HYBRID TECHNIQUE 
Governing Equations 

The steady-state thermal response of a structure or solid can be described by 
the following differential equation and boundary conditions; 

S£(T) = 0 in Q (1) 


and 


B(T) = 0 on r 


(2) 


where T is the temperature, Sf and 8 are differential operators, Q is the 
domain of the structure or solid, and T is its boundary. The differential operator 
includes nonlinear conduction, convection, and radiation terms. 

The application of the hybrid perturbation/Bubnov-Galerkin technique to the 
solution of equations (1) and (2) can be conveniently divided into the following two 
distinct steps: (1) generation of coordinate functions (or temperature modes) using 

the standard regular perturbation method and (2) computation of the amplitudes of the 
coordinate functions via the Bubnov-Galerkin technique. The procedure is described 
in detail subsequently. 


Generation of Coordinate Functions 

For the purpose of generating the required coordinate functions (or temperature 
modes), the governing differential equation (eq, (1)) is embedded in a single- or 
multiple-parameter family of equations of the form: 


5(T,q) = 0 


(3) 


for the single-parameter case and 


^(T,q^ ,q^) =0 (4) 
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for the two -parameter case, where q, , and q^ are normalizing or perturbation 
parameters. Extension to more than two parameters is straightforward and is not 
discussed herein. 

For nonlinear thermal problems, _it is convenient to choose the perturbation 
parameters such that the operators Q. and Q. can be decomposed as follows; 


5(T,q) = i^(T) + S£^(T,q) 


(5) 


and 


^(T,q^ ,q_^) 


= Sf (T) + (T,q, 

oo 11 




( 6 > 


where Sf and are JU^nOXVi differential operators which are independent of 

o oo 

q, q-| , and <^2f SP and are nonJLinZCUi operators in T and its spatial 

derivatives • 

The temperature T is represented by the regular perturbation expansion: 


n-1 

T = X T. (7) 

i=0 ^ 


for the single parameter case, and 


"i ^ . 

T = E Z qr" 


i 



T. . . 

i,:-i 


( 8 ) 


for the two-parameter case, where and are perturbation functions which 

represent temperature modes; and n is the total number of terms in the expansion. 
For the two-parameter case n = + 1 ) (n^ +2). The differential equations used 

in generating (or j) are obtained by substituting the temperature expansion 

(eq. (7) or eq. (8)), into the governing differential equation (eq. (3) or eq. (4)) 
and setting the coefficients of q^ (or q^q^ ) successively equal to zero. This 
leads to the following recursive set of equations (for a detailed description of the 
regular perturbation technique, see any of the monographs on the subject, e.g., 
refs. 7 and 8): 

i (T ) = 0 (9) 

o 0 

2'.(T. ) = R. (j > 0) (10) 

J 3 3 
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for the single-parameter case, and 


Sf (T ) 

oo 0,0 


0 


( 1 1 ) 


Sfl . (T. . ) = R. . (i > 0, j > 0 or i > 0, j > 0) (12) 

11 1/1 ID 

for the two-parameter case where Sf ! , and 2(!! . are UjKKX/i differential 

^ o ' oo 3 13 

^ /V 

operators; R. depends on Tq, with A < j; and R. . depends on T. , with 
3 ^ ■*'1 X, K 

X<i/ k<j or X<i, k<j. 

Note that whereas the original differential equation (eq. (3) or eq. (4)) is 
nonlinear, the recursive set of equations (eqs. (9) and (10) or eqs* (11) and (12)) 
is aUi ZJjitoji. 

For prescribed nonzero values for the boundary temperatures or their spatial 
derivatives, the zeroth-order perturbation function Tq (or Tq q) is used to sat- 
isfy these nonhomogene ous boundary conditions. All the higher-order perturbation 
functions, solutions of equation (10) (or eq. (12)), satisfy only homogeneous 
boundary conditions. 


Computation of Amplitudes of Coordinate Functions 


The perturbation functions (or Tj^^j) now chosen as coordinate func- 

tions and the temperature is expressed as a linear combination of these functions as 
follows: 


T 


n-1 

E 


i=0 


6. T. 

1 


for the single-parameter case, and 




T = y' T] . T. . . 


(13) 


(14) 


for the two-parameter case where and j are unknown parameters which repre- 
sent amplitudes of the coordinate functions (or temperature modes) and ^ 

and n equals the total number of modes. In the two-parameter case. 


n 



+ 1 ) (n^ 


+ 2 ) 


The parameters (or ^i,j) obtained by applying the Bubnov-Galerkin 

technique to the governing differential equation (eq. (3) or eq. (4)). The resulting 
set of n nonlinear equations can be cast in the following form: 



4-iTi,q 


dQ = 0 


(15) 


6 


or 



(b. . .T. . . 

x,3-i 



T 




dQ 


0 


(16) 


Note that for prescribed nonzero boundary conditions, the parameter cpQ (or (\>q q) 
is used to satisfy these boundary conditions; therefore, the number of free parame- 
ters reduces to n-1 • 


Comments on Selection of Coordinate Functions 

The chosen set of coordinate functions has the following two properties: 

1. They are linearly independent and span the space of solutions in the neigh- 

borhood of the point of their generation. Therefore, they fully character- 
ize the nonlinear solution in that neighborhood, 

2, Their generation, using the regular perturbation technique, requires the 

solution of a recursive set of t^ntCUl differential equations. 

The first property is necessary for the convergence of the Bubnov-Galerkin 
approximation, Ihe second property enhances the effectiveness of the proposed hybrid 
technique for solving nonlinear thermal problems. 

The mathematical complexity of the governing differential equations used in 
generating the coordinate functions builds up rapidly with the increase in the number 
of these functions. Therefore, it is convenient to generate only a few coordinate 
functions and augment them, for example, by the first eigenfunction of the zeroth- 
order perturbation equation. The performance of the augmented set of coordinate 
functions will be discussed in the next section. 


NUMERICAL STUDIES 

To evaluate the effectiveness of the proposed hybrid analysis technique, three 
nonlinear steady-state thermal problems were solved. For each problem, the solutions 
obtained by the hybrid technique were compared with exact solutions (whenever avail- 
able), the regular perturbation expansions, and other numerical approximations. 

Herein the results of three typical steady-state thermal problems are discussed. The 
three problems are (1) two-dimensional conduction in a square plate, (2) a one- 
dimensional conduc ting-con vec ting fin with variable heat-transfer coefficient, and 
(3) a one-dimensional conducting-convecting-radiating fin. In the first and third 
problems, the thermal conductivity is assumed to be temperature dependent. The three 
problems were analyzed by using the regular perturbation technique in reference 9. A 
summary of the differential equations and boundary conditions for the three problems 
is given in table I. 

To assess the accuracy of the solutions obtained by the different techniques for 
each of the three problems, the temperature is evaluated at a large number of points 
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in t±ie solution domain. llien a vector of temperature errors {AT} is computed as 
follows: 


{AT} = {T> 

e 


- (t) 

approx 


where {t}^ and are the vectors of QXOict and approximate temperatures, 

respectively. 


A weighted Euclidean norm of {AT} is used as a global error measure, namely 


{at}^ (at) 

{T}^ {T> 
e 


where superscript t denotes transposition, and m is the number of points at which 
the numerical values of the temperature are evaluated. For the three problems ana- 
lyzed, m was chosen to be 49, 30, and 30, respectively. 


Two-Dimensional Conduction in Square Plate 

The first problem considered is that of steady-state thermal conduction in a 
thin square plate with prescribed boundary temperatures. (See fig. 1.) The thermal 
conductivity is assumed to vary linearly with temperature. An exact analytic solu- 
tion based on the Kirchhoff transformation and a perturbation solution were presented 
in reference 9. As in reference 9, the perturbation parameter q was chosen in the 
present study to be the nonlinear conductivity coefficient y. (See table' I.) 

Two sets of coordinate functions are considered in the hybrid technique. The 
first set uses the two terms of the perturbation expansion given in reference 9. In 
the second set, the two coordinate functions of the first set are augmented by the 
lowest eigenmode of the zeroth -order perturbation equation. The explicit forms of 
the coordinate functions are listed in the appendix. Since one coordinate function 
(zeroth-order perturbation solution) is used in satisfying the nonzero prescribed 
boundary condition, the first set of coordinate functions allows the use of a single 
free parameter and the second set has two free parameters. The numerical values of 
the amplitudes of the coordinate functions in the hybrid technique are given in 
table II for the two cases q = 1 and q = 3. Note that the amplitudes of the first 
two-coordinate functions in the regular perturbation technique are equal to 1 and q, 
respectively. 

An indication of the accuracy of the solutions obtained with the hybrid tech- 
nique and the perturbation method is given in figure 1 . The error norm for four 
different values of the nonlinear conductivity coefficient (perturbation parameter), 
namely, q = 0.5, 1.0, 3.0, and 5.0, is shown in figure 2. An examination of fig- 
ures 1 and 2 reveals the following: 

1 . As expected, the accuracy of the perturbation solutions deteriorates rapidly 
with the increase in the perturbation parameter q (fig. 2), particularly for 
q > 1.0. For q = 3.0, the perturbation solution was considerably in error 
(fig. 1), and for q = 5.0, the error norm for the perturbation expansion could not 
be shown in figure 2. 
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2. Ilie solutions obtained by using the hybrid technique with two-coordinate 
functions (one free parameter) were reasonably accurate for q < 3.O. For q > 3,0, 
the hybrid technique predicts qualitatively the correct temperature distribution. 

3. The accuracy of the solutions obtained by the hybrid technique was considera- 
bly improved when the two-coordinate functions were augmented by the eigenmode. This 
is particularly true for q > 3.0. (See figs. 1 and 2.) 


One-Dimensional Conducting-Convecting Fin With 
Variable Heat-Transfer Coefficient 

The second problem is a straight conducting-convecting fin of length L, cross- 
sectional area and perimeter c, exposed on both sides to a free convective envi- 

ronment of temperature T , (See fig. 3.) The boundary conditions are a constant 
base temperature and an adiabatic tip. The thermal conductivity \l is assumed to be 
independent of the temperature. The convective heat- transfer coefficient is taken to 
be of the form: 


h = ( 18 ) 

b 

T - T 

where 0 = — — is a normalized temperature defined in terms of the ambient tem- 

b ^a 

perature T^ and the fin base temperature Tj^; p is a small parameter 

(p = 0.25 and 0.33 for laminar and turbulent conditions, respectively (ref. 9)); 

and hj^ = where y is a constant. 

Two different choices were made for the perturbation parameter q. The first 
choice is the same as that of reference 9, namely q = p. The second choice is 


q 


h cL^ 
b 

feA 


2 

where X is a convection-conduction fin parameter (table I) , The latter choice 
resulted in significantly simplified expressions of the coordinate functions. For 
case 1, q = Pf two-coordinate functions were used, and for case 2, q = X^, four- 
coordinate functions were generated. This corresponds to a single free parameter in 
case 1 and three free parameters in case 2. The expressions of the coordinate func- 
tions for the two cases are given in the appendix. 

7\n indication of the accuracy of the solutions obtained by the hybrid technique 
and the regular perturbation method are given in figure 3 for X = 1.0 and 2.0 and 
p = 0.33 and 1.0. The latter value of p has no physical significance and was 
selected in order to amplify the effect of the magnitude of the perturbation parame- 
ter on the quality of the solutions. The standard for comparison was taken to be the 
finite-element solution obtained by using a uniform grid of 15 three-noded finite 
elements with quadratic Lagrangian interpolation functions for the temperature. As 
can be seen from figure 3, the accuracy of the perturbation solution is very sensi- 
tive to both the choice and the magnitude of the perturbation parameter. For 
q = P, the two-term perturbation expansion is accurate for p = < 0.33 but becomes 
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quite inaccurate for p = 1.0. On the other hand, for q = , the four-term 

perturbation expansion is grossly in error for all \ > 1.0. The perturbation 
solutions for X = 2 could not be shown in figure 3. By contrast, the accuracy of 
the solutions obtained by the hybrid technique were found to be insensitive to the 
choice of the perturbation parameter. The solutions obtained by using q = p 
and q = were equally accurate and were almost indistinguishable from the finite- 

element solutions. 


One-Dimensional Conducting-Convecting-Radiating Fin 

The last problem considered is that of one-dimensional conduction in a straight 
fin of length L, cross-sectional area A, and perimeter c. Heat transfer from the 
surface of the fin involves both convection and radiation. The thermal conductivity 
is assumed to vary linearly with temperature and the convective heat-transfer coeffi- 
cient and the surface emissivity are held constant. The boundary conditions con- 
sisted of a constant base temperature and an adiabatic tip. 

The solution to this problem using a two-parameter regular perturbation tech- 
nique was presented in reference 9. The two perturbation parameters were chosen to 
be the nonlinear conduction coefficient radiation-conduction parameter ji. 

(See table I.) The six terms of the perturbation expansion of reference 9 were used 
as coordinate functions, and since one coordinate function was used in satisfying the 
prescribed nonzero boundary condition, only five parameters are left as free parame- 
ters. The expressions of the coordinate functions are given in the appendix. 

The accuracy of the solutions obtained with the hybrid technique and the two- 
parameter perturbation method is shown in figure 4 for three different values of the 
nonlinear conduction coefficient and the radiation-conduction parameter (perturbation 
parameters q-^ and ^respectively, q-j = = 0.5, 1.0, and 3.0). The standard 

of comparison is taken to be the finite-element solution using a uniform grid of 
15 three-noded finite elements with quadratic interpolation functions for the temper- 
ature. As can be seen from figure 4, the solutions obtained with the perturbation 
method for q^ = ^ grossly in error. By contrast, the solutions obtained 

with the hybrid technique are in close agreement with the finite-element solution 
except for the small temperatures near the tip. The convergence of the solutions 
obtained with the hybrid technique with the increase in the number of coordinate 
functions is contrasted with the convergence of the perturbation technique in fig- 
ure 5. The error norms for the six- term perturbation expansion with q^ = q^ = 3.0 
and 5.0 could not be shown in figure 5. As can be seen from figure 6, the addition 
of the eigenmode to the coordinate functions resulted in reducing the errors near the 
tip. For X= 5.0, and q^ = <52 ” 3.0, the maximum error in the small temperature at 
the tip obtained by using six-coordinate functions was 18.9 percent. This error 
reduced to 10.6 percent when the eigenmode was added. 


POTENTIAL OF THE PROPOSED HYBRID TECHNIQUE 

The proposed hybrid analysis technique appears to have high potential for solu- 
tion of nonlinear steady-state thermal problems. The numerical studies conducted 
clearly demonstrated the accuracy and effectiveness of the technique. In particular, 
the following two points are worth mentioning: 
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1 • rrhe proposed hybrid technique can be thought of as either of the following: 

(a) A generalized perturbation method in which (1) the perturbation expansion 
of the temperature contains free parameters, rather than fixed coefficients and 

(2) the perturbation parameters need not be small. Since the accuracy of the solu- 
tions obtained with the hybrid technique is insensitive to the choice of perturbation 
parameters, they may be introduced artificially to simplify the form of the recursive 
set of differential equations used in evaluating the various-order perturbation 
solutions • 

(b) An extended Bubnov-Galerkin approach with the coordinate functions gener- 
ated by using the standard regular perturbation technique, rather than chosen 

a priori. 

2. Ihe hybrid technique presented herein is the analytic counterpart of the 
reduction method presented in reference 6. In the reduction method, the initial 
discretization is done via finite elements, then the vector of fundamental unknowns 
is expressed as a linear combination of a small number of global-temperature modes or 
basis vectors and the Bubnov-Galerkin technique is used to compute the coefficients 
in the linear combination. The primary objective of using the reduction method is to 
reduce considerably the number of degrees of freedom in the initial discretization, 
and hence, reduce the computational effort involved in the solution of the problem. 

By contrast, the objectives of the foregoing hybrid technique are (a) to extend the 
range of validity of the regular perturbation method by removing the restriction of a 
small perturbation parameter and (b) to enhance the effectiveness of the Bubnov- 
Galerkin technique by removing the arbitrariness in the selection of the coordinate 
functions . 

3. The hybrid technique can be applied, in conjunction with other numerical 
discretization techniques, to the nonlinear thermal analysis of practical structures 
with complicated geometries. To accomplish this, a partitioning scheme can be used 
with the structure divided into substructures. The substructures with simple geome- 
try are analyzed by using the foregoing hybrid analytical technique, and the sub- 
structures with complicated geometry (e.g., near the boundaries or interfaces) are 
analyzed by using numerical discretization procedures. 


CONCLUDING REMARKS 

A hybrid analysis technique based on the combined use of regular perturbation 
expansion and the classical Bubnov-Galerkin approximation is presented for predicting 
the nonlinear steady-state temperature distributions in structures and solids. The 
application of the technique to the solution of nonlinear thermal problems can be 
conveniently divided into the following two stages: generation of coordinate func- 

tions (or temperature modes) using the standard regular perturbation method and 
approximating the temperature by a linear combination of these modes. The classical 
Bubnov-Galerkin technique is then used to compute the coefficients of the linear 
combination (amplitudes of the temperature modes). 

Three numerical examples demonstrate the effectiveness of the hybrid technique 
for the solution of nonlinear steady-state thermal problems. The three problems are 
(1) two-dimensional conduction in a square plate, (2) a one-dimensional conducting- 
convecting fin with variable heat-transfer coefficients, and (3) a one-dimensional 
conducting-convecting-radiating fin. In the first and third problems, the thermal 
conductivity is assumed to vary linearly with temperature. 
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The results suggest several conclusions relative to the choice of the coordinate 
functions and to the effectiveness of using the proposed hybrid technique in nonlin- 
ear steady-state thermal problems. These conclusions are as follows: 

1 . The proposed hybrid technique exploits the best elements of the regular per- 
turbation method and the Bubnov-Galerkin technique as follows: 

(a) The regular perturbation method is used as a systematic and general 

approach for generating coordinate functions. 

(b) The Bubnov-Galerkin technique is used as an efficient procedure for mini- 

mizing and distributing the error, in the temperature approximation, 

throughout the domain. 

2. The proposed hybrid technique extends the range of applicability of the per- 
turbation method and enhances the effectiveness of the Bubnov-Galerkin technique. It 
also alleviates the following major drawbacks of the classical techniques: 

(a) The requirement of using a small parameter in the regular perturbation 

expansion is avoided. 

(b) The method provides a systematic selection of the coordinate function (or 

temperature modes) needed in the Bubnov-Galerkin technique. 

3. For a given number of coordinate functions, the accuracy of the solutions 
obtained by the hybrid technique can be improved by the addition of the first 
eigenmode of the zeroth-order perturbation equation to the coordinate functions. 

This is particularly true when the number of coordinate functions is small. 

4. The accuracy of the solutions obtained by the hybrid technique is insensitive 
to the choice of the perturbation parameter (s ) • Therefore, the parameter(s) may be 
introduced artificially to simplify the form of the recursive set of differential 
equations used in evaluating the various -order perturbation solutions (viz, the coor- 
dinate functions) . 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
March 21, 1983 


APPENDIX 


COORDINATE FUNCTIONS USED IN NUMERICAL STUDIES 

The expressions for the coordinate functions used in the three problems are 
given in this appendix. These expressions were generated (and checked) by using the 
computerized symbolic manipulation system MACSYMA (ref. 10). The different symbols 
used in the expressions are defined in table I. 


Two-Dimensional Steady Conduction in Square Plate 
The two coordinate functions are given by 


0^ = sinh C^/sinh n 

1 9 

" ^ 2 2 ** 1 

0 = ^ 2[m(m - 4)ii sinh %] [ (1 - cosh 2u)sinh m? /sinh mTi 

^ m=1 ,3, . . . ^ 


- (1 - cosh 2F^^)]sin m^^ 


The first eigenfunction for the zeroth-order perturbation equation is given by 


6e = 


Steady-State Analysis of One-Dimensional Conducting-Convecting Fin 
With Variable Heat-Transfer Coefficient 

o 

The coordinate functions for the two cases q = fB and q == \ are given by the 
following equations: 

Ca se 1 ; q - g 


0Q = sech \ cosh 


0 


1 


C^ cosh Xl + C^l sinh Xl - C3|a4 + 42 


1 3 



APPENDIX 


where 


1 2 / ,2 41 4 37 6 

= — sech \ 84 + 42\ + — ^ X + X + 

12 V I ^ JDU 

- \ tanh X sech \ In sech X 




= 


X sech X In sech X 


C 3 = - sech X 


F, = 1 - 5 

Case 2; g = \^ 

^ 

®2 " 

83 = 1(1 + p)f- ^(2 + 3p)C + |(1 + p)C^ + 2 p^'‘ - |(1 + 4p)C^ + j ^(1 + 4p)c' 


Steady-State Analysis of One-Dimensional Conducting- 
Convecting-Radiating Fin 

The six coordinate functions in this case are given by the following equations: 


0Q Q = sech X cosh Xt 

12 - - 
0 ^ = ^ sech \[cosh 2X sech X cosh - cosh 2 \^] 

u f 1 s 
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- / 4 

-1 + sech \ cosh X^(1 - -g 

5 


cosh 2X - cosh 4X1 


4 — 1 

+ ~ cosh 2XC + cosh 4X^ 


® 0,2 " i 


cosh sech^X cosh^2X - ^ sech X cosh 3X - ^ X tanh X^ 


'.(± 

^\3 


)*i 


1 


- cosh 2\5{— sech X cosh 2Xj + ^ cosh 3XC + ^ X? sinh XC 


■ e 


1 ,1 


{- 


sech^X(-cosh 2X sech X 

2X^ 


+ sech X cosh X^hr sech X cosh 2X 


(-1 


cosh 2X 


15 


cosh 4X 


1 03 13 9 . 

+ cosh 3X + cosh 5X + X sinh X 


'] 


- 1 sech X cosh 2X|^1 - cosh 2X - ^ cosh 4X^ 

- 3^ cosh 3X^ + sech X cosh 2X cosh 4X1 


cosh 5X^ " "Ij sinh XC 


©2 0 ” sech^X^ ^^^ sinh X^{sech X[sinh 3X1(9000 - 4000 cosh 2X 

-200 cosh 4X) + sinh 5X1(1080 - 480 cosh 2X - 24 cosh 4X) 1 


- 12 960 sinh 2X1 *“ 120 sinh 4X^ + 480 sinh 6X1 + 15 sinh 8X? 


22 680X^} -* cosh Xl(2l6 cosh 2X^ + 240 cosh 4X1 + 98 cosh 6X1 

6 


8640 

+ 3 cosh 8X1) + sech X 


sinh^X^ - cosh^X?)(45 - 20 cosh 2X 


- cosh 4X) 


1 5 - \ 

+ ~ cosh X5 I + C cosh X^ 
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where 

1 8 

C = T sech \{sinh \[22 680X + 12 960 sinh 2X + 1 20 sinh 4X 

57 eoox"^ 

- 480 sinh 6\ - 1 5 sinh 8X - sech \ sinh 3X(9000 - 4000 cosh 2X 

- 200 cosh 4X) - sech X sinh 5X(1080 - 480 cosh 2X - 24 cosh 4X) ] 

+ cosh X(1080 cosh 2X + 1200 cosh 4X + 440 cosh 6X + 1 5 cosh 8X) 

2 

+ sech XEsinh X(-54 000 + 24 000 cosh 2X + 1200 cosh 4X) 

+ cosh^X(17 280 - 7680 cosh 2X - 384 cosh 4X) ] - 21 600 cosh^X} 


and 


1 = 1 -^ 


The first eigenfunction of the zeroth-order perturbation equation is 


0 = COS 

e 2 


Note that the expression for 9 is different from that presented in reference 9, 

2/0 
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TABLE I, - DIFFERENTIAL OPERATORS AND BOUNDARY CONDITIONS USED IN NUMERICAL STUDIES 





TABLE II.- AMPLITUDES OF COORDINATE FUNCTIONS IN HYBRID TECHNIQUE 


Square plate with temperature -dependent 
conductivity; q = y 


thermal 


1 


Amplitude 

q = 

1 

q = 

3 

II 

c 

n = 3 

n = 2 

n = 3 

4'0 

1 .0 

1 .0 

1 .0 

1 .0 

4^1 

0.6856 

0.6033 

1 .2278 

0.9283 

• 

• 

• 

• 

• 


0.01 56 


0.0581 


*Amplitude of eigenfunction 0^. 





4.0 


] q = 1.0 



0 2-coordinate 2-coordinate 

functions functions and 


PERTURBATION eigenmode 

METHOD HYBRID TECHNIQUE 

Figure 2.- Accuracy of solutions obtained by perturbation method and hybrid 
technique. Square plate with temperature -dependent thermal conductivity 
(problem 1 ) • 





Figure 3.- Comparison of solutions obtained by perturbation method and hybrid 
technique for one-dimensional conducting-convecting fin with variable heat 
transfer coefficient. 




0 .2 .4 .6 .8 1.0 

hi'- 


Figure 4,- Comparison of solutions obtained by perturbation method and hybrid technique for 
one-dimensional conducting-convecting-radiating fin with variable thermal conductivity. 

X = 2.0; is nonlinear conductivity coefficient; q 2 is radiation-conduction fin 

parameter. 
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Figure 5.- Accuracy and convergence of solutions obtained by perturbation method and hybrid technique. 
One -dimensional conducting-convecting-radiating fin with variable thermal conductivity; q.| is 
nonlinear conductivity coefficient; P 2 is radiation-conduction fin parameter. 



Finite-element solution 
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